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We present a detailed analytical analysis of the propagation of Bogoliubov phonons on top of 
Bose-Einstein condensates with spatial and temporal step-like discontinuities in the speed of sound 
in the hydrodynamic limit. We focus on some features in the correlations patterns, in particular of 
density-density correlations. The application to the study of the Hawking signal in the formation 
of acoustic black hole-like configurations is also discussed. 
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>>i I. INTRODUCTION 



The properties of quantum systems in which some parameters experience instantaneous changes in time (or quantum 
quench, see for instance [lj]-|3l]), have been considered with applications in many areas of physics. 

Such step-like discontinuities (in time, and even in space) have been analysed, in particular, for what concerns the 
propagation of Bogoliubov phonons on top of Bose-Einstein condensates [3]-]^ in connection with features of acoustic 
black holes and cosmology using the so called gravitational analogy [loj . 

Within the same context, our aim is to study correlations (in particular density-density correlations), by developing 
r : an approximation scheme to analyse, from an analytic point of view, the Hawking signal in the formation of acoustic 
^j-j' black- hole configurations, as considered in [ll|, [l^ (see also (T3|). 

In this analysis, restricted to the hydrodynamical approximation, we present a detailed discussion of phonons 
propagation in backgrounds with step-like discontinuities in the speed of sound, spatial (section 4.1) and temporal 
' (section 4.2) together with some features in the correlations patterns. Then, in section 4.3 we shall combine the two 
previous cases and repeat the analysis by considering the formation of an inhomogeneous (step-like) configuration 
'■ out of a homogeneous one. The application to study the Hawking signal from incipient acoustic black holes in the 
QQ . hydrodynamical limit is considered in the final section 5, where it is shown that a structure similar to that found in 
can indeed be recovered . 

§ ■ II. THE BASIC EQUATIONS 



X 



In the dilute gas approximation, a Bose gas can be described |14i] by a field operator ^' with equal-time commutator 

[i>{t,x),i'\t,x')] = 6^x-x') (1) 

and time evolution given by 

ihBt'i' = f-^^V^ + Ve.t + ff*^^) * , (2) 

where m is the mass of the atoms, Vext the external potential and g the nonlinear atom-atom interaction constant. 
By considering linear fluctuations around the classical macroscopic condensate in the form 

*^*o(l + 0), (3) 
the condensate wavefunction ^I^o satisfies the Gross-Pitaevski equation 

ihdt'i'o = (-^^^ + Vext + gn) ^-0 , (4) 
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with n = 1^*01 the number density, while the Bogohubov-de Gennes equation for the fluctuations takes the form 



indt^^[-—W' --^VU + mc2(c^ + 0t)^ (5) 

\ 2m m Wq / 



the speed of sound being c — \f^- 



III. OUR MODEL 



We shah consider condensates of constant density n and velocity v (for simplicity along one dimension, say x). As 
explained in [T^ , nontrivial configurations are still possible provided one varies the coupling constant g (and therefore 
c) and the external potential but keep the sum gn + V^xt constant. In this way the plane-wave = y/ne**^"^"™"*, 
where v = is the condensate velocity and hw^ = ^ — ^ gn + Vext, is a solution of ^ everywhere. The transverse 
degrees of freedom are assumed to be frozen. 

To consider the hydrodynamic limit it is useful to introduce the two hermitean operators (density fiuctuation) 
and 9^ (phase fluctuation) via 



Inserting in ^ and expanding as 



aj7i](i, x) + a]r7,](i, x)* 
aje]{t,x)+a]e]{t,x)* 



we get the mode equations 



{dt + vdx)n^ + ^dlB^ = , 



(7) 



(8) 



The equal time commutator ([T]) becomes, in this representation, 

\n'(t,x),9^(t,x')\^iM(x- x') , (9) 
which after integration gives the modes normalization 



/ 



dx{n]B^^ - n^QW = i&jj> . (10) 



IV. HYDRODYNAMICS 

In the hydrodynamical approximation one probes distances much larger than the healing length ^ = . This limit 
is easily achieved in the density-phase representation ([5]), where in the second of eqs. we consider 

^(n^ - ^-3ln-) ^ , (11) 

n 4 n 

leading to 

{dt + vdx)n^ + ^dle^^Q , 

{dt+vdx)e^ + ^n^ . ^ ' 

Eqs. (|12p can be rewritten as a single second order equation for 9^ provided one algebraically extracts n} from the 
second 



n 



n 



idt + vdx)e' (13) 
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and substitutes it in the first to get the second order differential equation 



1 



[dt + vd^)^{dt + vd.,)e^ + dle^ = o 



(14) 



In the language of the gravitational analogy [T3l , 0^ becomes a massless and minimally coupled scalar field propagating 
in the curved acoustic geometry 



9i^y = — 
mc 



-{c^ -v^) -w 
10 
10 

1 



(15) 



and satisfying the Klein-Gordon equation 



(16) 



which corresponds (we remind that in our model n and v are constant) to eq. (jl4l) provided 9^ does not depend upon 
the transverse dimensions. The scalar product is given by, see eqs (jlOp and (1131) . 



= ^33' 



(17) 



We shall consider the simplest possible nontrivial configurations, namely step- like discontinuities of the speed of sound 
c in both x and t. In order to explicitly solve these problems we shall rewrite eq. ()14p in the more appropriate form 



dt 



dtO^ V 



a. 



(1 - ^)a,0i - ^dto^ 

c c 



= 



(18) 



Such configurations allow us to extract analytical results which can be meaningfully compared with the numerical 
analysis of more complicated regular configurations. 

To clarify the presence of the inside the square brackets in ([T5|. from we note that our four dimensional 
acoustic metric (1151) has J—q ~ - and that the inverse metric is 



-1 


— V 








— V 


— V 




























c2 



(19) 



i.e., the relevant t,r components carry an additional factor of -. 



A. Step-like discontinuity in x (stationary case) 



Let us consider a surface (say, a; = 0) separating two semi-infinite homogeneous condensates with different sound 
speeds, i.e. we consider c{x) = ciO{—x) + CrO{x). We shall consider v < 0, that is the fluid flows from right to left. 

To solve eq. in all space we need flrst to solve it separately in the left {x < 0) and right {x > 0) regions, then 
impose the appropriate matching conditions at the discontinuity. Looking at (|18p these are 

= , (20) 

where [ ] means variation across the jump at a; = 0. 

Eq. is easily solved in left and right regions in terms of plane-waves iJe"*"*^*'^^ with the standard linear 

dispersion relation {w — vk)'^ = (ck)'^. Being the model stationary the matching conditions (j20p are solved at flxed w. 
This motivates the use of the w expansion for the operator 9^ 

r+ca 

9\t,x)^ dw[alulit,x)+alulit,x)+h.c.] (21) 

Jo 



(1 - -^)d,9' - -^dt9' 
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in terms of right-moving 



and left-moving waves [22 



'iwt-^ik-u.x r, 



W 



V -\- c 



(22) 



.V r> —iwt+ik^x h. 



V — C 



These modes are ortho-normalizcd according to 



i I dx [ul,{dt -I- vdx)u^, - u^,{dt + vdx)ul,] = S{w - w') , 



and for infinite homogeneous condensates we have 



Ann 



dk 

d'W 



Aimw ' 



(23) 



(24) 



(25) 



where = w — kv is the frequency measured in the atoms frame. 

Our derivation of is based on the spatial integration of the commutator which is time independent. To 
further check that this scalar product does not depend on time we consider 



(26) 



(27) 

Integration by parts in the third and fourth terms and use of the boundary conditions (j20p for the surface terms shows 
that this quantity indeed vanishes. 

In the 1 region (a; < 0) the generic solution to (|14p at fixed w is thus written in the form 



Using the equation of motion (jlSp , the above expression can be rewritten as 



dt \ dx 



A-Kwu 



I — — — X a! , I — r — X a! 

e ^""i At + e At, 



whereas in the r (x > 0) region 



(28) 



Attwu 

We shall now impose the matching conditions pO| at a; = 0. They read 



(29) 



-i{Ai + Al) ^ ^r{Al + Al) 



= ^^{[(1-^)^ + 1 
They can be rewritten in matrix form 



A. 



(1 



or more explicitly 



(30) 



(31) 



1 1 
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A,, 



1 1 

1 1 



A^. 



(32) 
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By inverting Wi we get 



A., 



M, 



the scattering matrix Mgcat = ^Wr being 



Mscat 



1 f Cr + Ci Cr - Cl 



V Cr " C; Cr + Cl 



(33) 



(34) 



1. 'in' and 'out' basis and two-point function 

Let us construct the 'in' modes , see Figure [TJ 







FIG. 1: 'in' and 'out' modes. 



These are initially left moving modes from the right (m^'™) and initially right-moving ones from the left (u^^*"). 
Modes 

Consider a left-moving (v) wave from the right (r) region with unit amplitude, i.e. = 1. The scattering at 
X = will produce a reflected wave (right- moving u wave in the r region, with amplitude A^^) plus a transmitted one 
(left-moving in the left(l) region, Al^). Thus we have to solve the following system of equations 



Its solutions are 



satisfying the unitarity condition 







1 



Cr +Ci Cr - Cl 



2v^CVQ \Cr - Cl Cr + Cl 



Cl n ./ 2y/CrCl 



1 



R,Al^ ^ "' = T 

Cl Cr+ Cl 



(35) 



(36) 



(37) 
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We note that this construction is valid only when the flow is everywhere subsonic. Indeed, if the left region were 
supersonic, \v\ > ci, we would have two transmitted components instead of one because there also the 'right-moving' 
one (Ajj) would be dragged to the left by the fluid (i.e. ci + v < 0). Mathematically the system does not admit a 
unique solution since we would have three amplitudes to determine with only two equations. From the physical point 
of view we would need a smoother configuration, an acoustic black hole, with a transition region from the right to 
the left regions of nonvanishing width Aax (see the final section 5). 

Modes 

A right-moving (u) wave from the 1 region with unit amplitude, i.e. = 1, gives rise to a reflected wave (left-moving 
V wave in the I region, with amplitude A[,) plus a transmitted one (right-moving in the r region, A^). 
The amplitudes are obtained by solving 



1 



1 



that gives 



Cr - Cl 
Cr + Cl 



Cr + Cl Cr — Cl 
Cr - Cl Cr + Cl 



-R,Al 



Cr + Cl 





Al 



T 



(38) 



(39) 



Again, this construction is for subsonic fiows. Were the left I region supersonic, the initial u mode could not propagate 
to the right, but would be dragged to the left by the fiuid. 

Similarly, one can construct the 'out' modes which are represented in Figure [TJ 



The relations between the 'in' and 'out' modes 



-Ru, 



v,out 



allow us to determine the relation between the a, a} operators in the 'in' and 'out' decompositions 



+ 00 



dw 



aj:'"(°"*)M:;'''"(°"*Hi, a;) a:^,'™(°"*)<'"(°"*) (t, x) /i.e. 



(40) 



(41) 



Ray + Ta]T 



Since (|42[) only involve annihilation operators the vacuum state is unique \in) — \out). 
We can then write down the symmetric two-point function 

1 



{in\{9^{t,x),e^{t' ,x')}\in) = - {in\Q^(t,x)9^{t' ,x') + B^{i! ,x')B^{t,x)\in) = 



+00 



dw 



2 Jo 

where, explicitly. 



ul;''''^°'''^ {t, {t' , x') + (t, a;)u^'"(°"*)* {t' , x') + c.c. 



= uy{t, x) + Ruy{t, x) + Tul;\t, x) , 
= uZ''it,x) - RuZ'\t,x)+TuY{t,x) , 

with left and right moving components in the I and r regions 

"-'^''"^^ ^ \f^^^^'^'^^^^' <'-(t,x)=<-'(t,x)[Q^c 

Ul;'{t,x) = V K^'it^x) ^ui{t,x)[ci Cr] 



(42) 



(43) 



(44) 
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Considering for instance points situated, respectively, in the 1 (2; < 0) and r {x' > 0) regions we have 

h 



{in\{9\t,x),e\t',x')}\in) ^ 



duj 



r*<''(i,xX"(i',x') + 



T<''(t,a;X"(t',a:') + {R*T ^ RT*)uZ^Ht, x)uy* , x') + c.c. 



hmCrCl 

27rn(cr + q) 



ln[(i - t' 



V + Cl V + Cr 



V - Cl 



(45) 



The non vanishing contributions are of uu and vv type, while the uv correlations, being proportional to {R*T — RT*), 
vanish. 

In section III. A of [l^ we have considered dispersion effects for spatial step-like discontinuities in subsonic flows. 
There it is shown that the ^ — > limit is smooth and that the hydrodynamic modes give the main contribution to 
the correlators. 

Finally, we stress that from the point of view of the gravitational analogy the nontrivial scattering structure here 
outlined depends crucially on the fact that the geometry (|15p is four-dimensional. Were g^i, two-dimensional, the field 
9^ would be conformal and there would be no scattering between right-moving and left-moving components. 



B. Step-like discontinuity in t (homogeneous case) 



The other case of interest is the step-like discontinuity in t (say, at i = 0), separating two infinite homogeneous 
condensates: c{t) = CinO{—t) + CoutS{t)- 

To solve for all times means we need to write down the general solutions in both 'in' {t < 0) and 'out' {t > 0) 
regions and then impose the appropriate matching conditions for the problem at t = 0, which are (see psp ) 



m - 0, 







(46) 



In both 'in' and 'out' regions the solutions are plane- waves e and being the problem homogeneous the 

boundary conditions pS)) are solved at fixed k. Therefore we shall use the k decomposition for 9^ 



9\t,x) 



+ 00 



dk 



akUk{t,x) +aluk{t,xy 



(47) 



where the positive-frequency modes 



Uk 



(48) 



are right-moving, Wu = {v + c)k, when fc > 0, and left-moving, Wy — {v — c)k, as fc < 0. Such modes are ortho- 
normalized according to 



{uk,Uk') = i / dx 2 [K{dt + vdx)uk' - Uk'{dt + vdx)u*k] = S{k - k') , 



(49) 



I.e. 



This scalar product too is time independent, i.e. 



47r|ri|n 47r|fc|n 



(50) 



iuk,Uk')t2 - {uk,Uk')ti = i- 



dt / dxdt 



u*k\{dt + vdx)uk' - Uk'\{dt + vdx)ul 



. 



(51) 



For tit2 > this is trivial, whereas for tit2 < this is guaranteed by the boundary conditions (|46p at the timelike 
discontinuity. See also 0. 
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Note that the general sohition at fixed k 



(52) 



contains a positive-frequency and a negative-frequency component: when fc > (fc < 0) it describes a positive 
frequency right-moving (left-moving) mode and a negative frequency left-moving (right-moving) one. 
In the 'in' region we write 



t1 _ ikx I '^C; 

'in ^ 



47r|fc|n 



A,, e -|- /i,, e 



whereas in the 'out' region 



ikx / ^^out 



A'n:\k\i 



Aout —ik{v—Cout)i I Aout ~ik{v+Cout)t 



(53) 



(54) 



The matching conditions (l46l) read 



r{[(^^-Co.t)-^^]^r 
i^{[(t^-c„)-«]A- 



[(«- 



(55) 



In matrix form 



or more explicitly 



A \ I A' 



Wont ( ^:^ut j=W,n[ ) (56) 



^out ^out 

-1 1 



j^OUt 
A out 



By inverting Wout we get 



A out 

J^out 



the Bogoliubov matrix Mboq = ^o^i^^^?^ being 



M 



Bog 



^inCout 



-1 1 



~ ^out C-i't 



Cout 
C-out 



J^ut 



(57) 



(58) 



(59) 



1. Connecting the 'in' and 'out' basis 



The 'in' and 'out' basis are easily identified. They are constructed in terms of 'in' and 'out' positive-frequency 
modes which are left-moving /c < 0, and right-moving u™''""*'', fc > 0. To connect them, as depicted in Figl2] 

for subsonic flows, we use the Bogoliubov matrix (I59p . Were the 'out' region supersonic, both the modes would move 
to the left. 

Let us construct explicitly the 'in' modes. 
Modes u^" (fc < 0) 

A left moving 'in' mode has A™ — 1, A™ = 0. The coefficients and A""* are found by solving 

(^'Q I [ Cin 4" ^out C-in Cout \ ( ^ \ (C(W 



FIG. 2: 'in' and 'out' modes. 



The solution is 



A out _ 
— 



^out 



out _ '-"I 



* J^OUt 



^out 



These coefficients satisfy the unitarity condition 

i^rT - i^rT ^ lap - = 1 , 

the — sign meaning that the are associated to negative frequency and negative norm modes. 



(61) 



(62) 



Modes u^" (fc > 0) 



A right-moving 'in' mode has A™ = 0, A™ = 1. From 



Cm ~t~ ^out ^in ^out \ ( ^ 



we get 



^out 



A out _ ^in + Cout _ ^* 
I — ^ ■ - — ^ 



2-^ CinCout 



(63) 



(64) 



2-^ CinCout 

The 'out' modes can be similarly constructed, see FigO 

The relation between the 'in' and 'out' modes can be written compactly 

and it allows us to determine the relation between the annihilation and creation operators in the 'in' and 'out' 
decompositions 



/ + 00 
dk 
-oo 



in[out) in{out 



a,„ 'u 



P""^'(t,a;) + /i.c. 



(65) 
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i.e. 

^ out ^"in o*^mt fa£i\ 

flfe = a ftfc - P a_fc • (66) 

We see that l|66p involves both annihilation and creation operators, so the two decompositions are inequivalent and 
in particular their vacuum states are different \in) ^ \out). Up to a normalization factor these two states are related 
by the formula 

\in) - e"^ \out) , (67) 

showing that the \in) state contains correlated pairs of 'out' particles with opposite momenta (i.e. one right-moving 
and the other left- moving). These features show up in the correlations. 

Again, we point out that in the language of the gravitational analogy all the nontrivial features here discussed 
are a consequence of the four-dimensional character of the acoustic metric, and that they would be absent if g^^ were 
two-dimensional and, consequently, the field 6i conformal. 



2. Two-point function and density-density correlations 



Zero Temperature 

To calculate {in\0^{t,x)0^{t' ,x')\in) in the 'out' region, let us expand 6i in the 'out' basis 

/ + 00 
dk [a°^'ul''\t,x)+h.c] , 

and use the Bogoliubov transformation This gives 

{in\{e^{t,x),9^{t\x')}\in) = 
I dk \ {a*ul''* - I3u°^^*) {t, x) {aul""^* - /3*u°^^) {f , x') + c.c. 



hmco 



47m 



2c,„c„„t 



In A^^'A"^ 



t A out I (cL+cLf) 



uv' u'v ' 2Ci„Co„t 



It^ a out A out 



where we have defined 



l^Z' - (^ + Cout)t + [cout - v)t' ~{x~ x'), A°y* = AZ*>{t, X ^ t', x') , 
A°:* ^{v + c,.,t){t t') -{x- x'), AZJ = (c - v,^,)it - t') + {x- x') . 

The quantity of interest in the experiments is the normalized, one-time density-density correlation function 

G^''Ht;x-x') = \{in\{h\t,x),n^{t',x')}\in)\t=t' , 
that using (|13p can be calculated by performing derivatives from the above results 

2 



G(2) = 



1 



{dtdt.^vdtd.,,^vdt.d.,+v^d^d^\ (^n|{0l(^,a:),0l(^',x')}|m)|t=t, . 



The computation is straightforward and we get the following contribution 



G(2 



2 2/ 

2c,„Co„t \^ 



(2c„„tt-(x-2;'))^ + (2c„„tt-(a;'-x))^ 



(68) 

(69) 
(70) 

(71) 

(72) 
(73) 

(74) 
(75) 



The interesting features are in the time dependent contributions, coming from the a/3 terms in (|69p §3], and whose 
physical interpretation follows that of Eq. ((67|) . At t = and everywhere in space correlated pairs of particles 
with opposite momentum are created out of the vacuum state, with velocities v — Cout (left-moving) and v -1- Co„t 
(right- moving). At time t such particles are separated by a distance \x — x'\ = 2co„tt, which is indeed the correlation 
displayed in (I75|) j24j. As emphasized in the fact that the corresponding peak is infinite is due to the sudden jump 
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approximation (step-like discontinuity) employed in this analysis. A way to smooth out the transition from the 'in' to 
the 'out' regions within a time interval at and produce a finite quantity which is in good agreement with the numerical 
results is to multiply the integrand of Eq. by e~ , which introduces a cutoff at k^ax ^ ^ • This means that 
the relevant physics is in the low k (< kmax) modes for which the sudden transition is a good approximation to the 
real (smooth) situation. 



Finite temperature 

We shall now consider the case in which the initial state is not the vacuum \in) but thermal, e.g., it is described by 
a thermal density matrix at temperature T 



+00 



(76) 



N=0 



Therefore we shall need to calculate 



(77) 



To perform this calculation, starting from the decomposition (|68p and the Bogoliubov transformation (j66p we need 
now to transform to the 'in' w basis 



/'+00 



dw 



where, we remind, the fc > and k < 'out' modes are 



* out _ a,, out* 

" ^k=uj/{ci„+v) P^k=-w/ic,„+v) ^u,zn 
\/v + Cin ™ 



'mCput ^-ik(v+Cout)t+ikx out _ 



"mCput ^-ik{v~c„ut)t+ikx 



(78) 



(79) 



-(1) -(1) 

and CLw = "^^^j '^w = ■ Inserting (1751) into ([77]) as an intermediate step we get the following expression 



. . h 1'+°° 



'ifli = - / dw{l + 2\l3T 



/o 



2) (F{t, x)F* {t', x') + G{t, x), G* {t',x') + c.c.) , 



where \(3t? = F ^ """^-^'-"^--^Tf"^"---^'-"^"' and 



G 



"'^"fc"-^./(c^„^^,)-/^"g = ^/(c^„-^,) 



Comparing with the zero-temperature case we see that the 
the existing correlations through the factor 1 -f 2|/3tP, and not to ere; 
from (|79p . and evaluating the integrals we obtain the final expression 



effect of the initial thermal population is to stimulate 
to create new ones. Writing out explicitly the modes. 



Trpth0\t,x)0\t',x') = - 



AtTU X 2CmCout 



sinh 



sinh 



nkgT A°J 



out 

v 



+ 



^2 1 2 
^in ' ^out 



sinh 



out 
vv' 



h Cin—vJ_ 



(80) 



(2) 

From this result the calculation of , using ([71)1 . gives 



;x — x 



g2 , „2 



p2 _ 2 

in ^out 

^^in^out 



iTTnmCout I 2Cm Cout 

1 



52 



.sinh2[A(x - x')] ^ sinh2[B(.T - x')] 
1 



sinh2[A(2co„ti - {x - x'))] ^ sinh2[yl(2co„ti + {x - x')) 

1 



sinh2[B(2c, 



^ 2 

■outt - {x - x' j)] sinh [B{2coutt + {x ~ x'))] 



(81) 
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where A = ^^^^'^^ and B = Tj^^fz^- Comparing with (|75|) . we see that indeed the peaks structure is the same, and 
in particular the time-dependent ones receive the following finite contribution 



^q{2) ^ TrjkBTf (c^„ - cg„t) {v^ + cf„) 

T,\x~x'\=2c,^tt l2nmCout 2cinCout (W + Qn)^(Cm - W)^ ' 



which in the more realistic case of a transition within a time at between the 'in' and the 'out' regions and v — has 
been shown to be a good approximation to the numerical results for low temperatures T < h/kscrt 9]. 
Finally, it is noted in section IV. A of 15] that dispersion effects do not change qualitatively the hydrodynamical 
features discussed here. 



C. Simulation of the formation of non-homogeneous backgrounds using step-like discontinuities 



We shall now combine the analysis of the two previous subsections to construct more complicated time-dependent 
backgrounds. Consider an initial infinite homogeneous condensate with sound speed Ci„ turning into a nonhomoge- 
neous configuration well described by two semiinfinite (left and right) homogeneous condensates with sound speeds, 
respectively, ci and Cr- Technically the model we use consists of a discontinuity in t followed by a discontinuity in a;, 
i.e. c{t, x) = CinO{—t) + 9{t){ci9{—x) + CrO(x)). For simplicity we shall consider the case in which Ci„ = Cr, represented 
in FigH 



C ; 



in 




FIG. 3: Temporal formation of a spatial step- like discontinuity (cr — dn). 



Such a construction is the closest, in spirit, to the time-dependent formation of an acoustic black hole as considered 
in [ll| . 12 1, the important difference being, of course, that our hydrodynamic analysis can only be made for subsonic 
fiows. In addition, a realistic realization of this scenario requires non vanishing transition regions in both t and x of 
width at and ax- 

To solve (|14p for all x and t, first we write down the general solutions in the three different regions , 'in' (t < 0) 
and 'out' left (i > 0, x < 0) and right (t > 0, x > 0), and then impose the matching conditions (|46p and ((20)) . The 
former are performed at fixed k, while the latter at fixed w. Therefore both k and uj decompositions (|47l) and (PT|) 
will be needed, the respective modes and operators being related by (see also [Tsj ) 



Uk{<0) = Vc~V , Mfc(>0) = VcTv 

afe(<o) = Vc-v , afc(>o) = \/c + v . 
To solve (|46)) at t = consider the general solutions of (fT4l) at fixed k in the 'in' region 



(83) 
(84) 



47r|A;|n 



-ik(v- 



in —ik(v-\-c.in)t 



(85) 
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and in the 'out' left and right regions 

6)1 (fc) = e'^'-' 



mci 
47r|fc|n 



47r|fc|n 

For X < the 'in' and I amplitudes are related through the Bogoliubov matrix 



(86) 
(87) 



1 



2yfc~ci \ Cin - Cl Cin + Ci 

At x > the matching is trivial because we chose Cm — Cr, therefore 



AT _ Atn AT _ Am 



(89) 



The matching equations (UHl) at a; = for alH > relate the general solutions at fixed w in the I and r regions 



via the scattering matrix 



A!, 



4TTwn 



e 



1 



Cl 



2Jc^l \ Cin - Cl Cin + Cl 



A^ 

Ar, 



(90) 
(91) 

(92) 



In the next subsection we shall construct the explicitly the 'in' modes and study their time evolution in the t > 
region. This analysis is also useful to clarify what happens at the point x = t = 0, belonging to both the spatial and 
the temporal step-like discontinuities, and to the boundary conditions satisfied as one approaches it. 



1. Relating the 'in' and 'out' basis 



The way the 'in' k basis is related to the 'in' and 'out' w basis {t > 0) is shown in Figs. |4] and O where the 
Bogoliubov and scattering coefficients are 

* Cin + C; Cin Q „ Q ^ ^WCinCl /nn\ 

a ^ „ ^ : -P = „ ^ — , R = — , T = . (93) 

^yCinCl L^CmCl Cin ~r Cl Cin ~r Cl 

These pictures are constructed by combining those of Figs. [T] and [2] for the two types of discontinuities. 
Left-moving 'in' modes 

The left moving 'in' modes produce at t = 0^ three different components, two due to the discontinuity on the x < 
side and the third, freely propagating at x > 0. Only two of them scatter at the a; = discontinuity, the other one 
being disconnected from it (it propagates outside the future sound cone of the origin). Explicitly, 

f^r~ ^-ik{v-Ci„)t+ikx r^(n( ^\^* ^-ik{v-ci)t+ikx p^-ik{v+ci)t+ikx 

V - Cl 

+0(t ^)[pJl*e~'Hv+c,}t+^k!i±^x ^ ^^^^k{v-c,^)t+^k^x.y ^ 

V - Cl 

V + Cin 

We see that in the t > region there are several ways to get to the point x ^ t — 0. If we approach it from outside its 
sound cone from the left (the right part being trivial) we have the usual components arising from the temporal step- 
like discontinuity satisfying the boundary conditions (1461) . If instead we approach it from inside the sound cone the 
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FIG. 4: Evolution of M|"(fc < 0) modes. 




FIG. 5: Evolution of ui"{k > 0) modes. 



existing components there satisfy the boundary conditions (|20|) . Note that the mode (I96|) is continuous at x = t = 0, 
irrespective of how we get to it, and along the sound cone generators t — -^fz:^-, t — as one can easily verify from 

the relations a* = 13 R* + T and i? - l3T* =0. 

In the next subsection we shall perform the analysis of the correlations patterns and look for stationary signals (like 
the Hawking signal, see [ll| ) . For this reason we shall focus our attention to the interior of the future sound cone of 
the origin where we have 

— p\/v + cj{—H u , , , + 1 u , , , .) + \/ Cr — v(i u , N, + Ku , s, . 95 

^ ' w — ~k[v+ci) w— — k{v+ci}^ I V f V w—{v~Cr}k w—{v—Cr)k' ^ > 

We mention that with this restriction we are neglecting transient effects such as those arising from the temporal 
step-like discontinuity (see ([75)) ). see the numerical results in fl^]. Such effects are absent also in |llj| because there 
we considered the conformal field approximation for 9^ . 



Right-moving 'in' modes u^">q) 

In this case of the three components produced at t 



— 0+ only the one proportional to a* propagates inside the 
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future light-cone of the origin and scatters at the spatial step-like discontinuity. The explicit time evolution is 



-ik(v+Ciri)t+ikx 



-i{e{- 



V - Cl 



P'^^-i^('"-(^l)'t+ikx _|_ ^* ^-ik(v+ci)t+ikx 



V - Cl 



' V + Crn 



- t)e 



-ik{v-\-Cin)t-\-ika 



}• 



(96) 



The relations — /3 — —a*R, a*T = 1 guarantee continuity at a; = t = and along the sound cone generators 
t = ^TT^, t = . The boundary conditions involving derivatives in p6)) and pop are satisfied at the origin by 

approaching it from, respectively, outside (from the left) and inside its sound cone. 
Restricting again to the interior of the sound cone of a; = t = we have 



■in , * / ; uAn * / ; / 7-) v,out , rr u^out \ 

k{>{)) V t w — {v+ci)k ^ w — [v+ci)k w — {v+ci)k^ 



(97) 



Let us now construct the explicit connection between the a, operators of 'in' k decomposition {t < 0) 



+00 



dk 



arul-{t,x)+a"^lu^l\{t,x) + h.c. 



where we have explicitly separated the positive k (left-moving) and the negative k (right-moving) contributions, and 
the 'in' and 'out' w ones in the t > region 



+00 



e\t,x)^ / dw aj;'"'(°"*)<"(°"*)(t,a;)-Ka^'™(°"*)u^'"(°"*)(i,a;)-H/i.c. 
Jo L 

By inserting the first lines of eqs. (|95|). ((97)) into (l98l) and changing from fc to w we have 



r+oo p(a*afe-/3*alj|fe=_™_ 
71 = / dw 
'0 



a~k= 



„.u,zn _^ <;™ + h.. 



from which 



^/V + Cl VCin - V 

{a*dk -/3*al^)|fc= » 



a~k=- 



^ci + V 



(99) 



(100) 



(101) 



To change from the 'in' to the 'out' w basis we can use the second lines of (^5]). (^7]) into ([M)) or, directly, and 
dTUTl) : 



;-U,OUt _ p 



a-k=- 



T 



{a*du~P*d'_i,)\k=^ 

\/Cin " V ' ^Cl + V 

R — 



V Cin — V 



VCl+V 



(102) 



2. Correlations 



To evaluate in the t > region the two-point correlation function of 0i in the 'in' vacuum we can use the 'out' w 
decomposition 



+00 



9\t, x)^ j dw [a^'°"*u^'°"*(i, x) + C™*u^'™*(t, x) + h.c 



(103) 
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and then the relations (|102[) . This gives 



{in\§\t,x)§\t\x')\in) = 



2 7o 



Alternatively, we can use the 'in' w decomposition 

f +00 



(104) 



(105) 



from which using (jlQip 



+00 



{in\e\t,x)e\t',x')\in) ^ - / du. u;'*"(x,0<'"*(2:', <') + (! + 2|/3n<'™(x,i)<'™*(2:',0 + c.c. . (106) 



Given the relations (|40| between 'in' and 'out' w modes, the two expressions (|104p and (|106p are clearly equivalent. 
With respect to the stationary case P5)) {(3 — 0), we see in (|106l) the different weight the 'in' v and u modes have, 
modifying the correlation picture. In particular, considering points respectively in the left (say, a; < 0) and the right 
{x' > 0) region we have (see (|43|)) 

fc poo 

{tn\§\t,x)§\t',x')\zn} - - / dco T*{l + 2\(3f)u:('{t,x)uY*{t',x')+ 



Tu 



:;-'(i,x)<"(i',x') + ( 



R*T - RT 



-{l + 2\p\^))u:j\t,x)uY*{t\x') + c.c. 



hm 
Aim 



Cl 



ln{t - t' 



V + Cl V + c 



■in »^2n 



Cl 



V - Cl V - a, 



^ Cl) 



V - Cl V + Cv 



where, comparing with (|45p. we see here the different contributions coming from the uu and vv sectors and the 
appearance of a new term, a nonvanishing uv correlator proportional to (i?T*)|/3p. 
We can also compute the left-right density-density function 



(jnCinClY 



{dtdv^vdtd^.^vdvd.,^v'dA\ (^n|{^l(^,x),^l(^',x')}|^n)|^=^' 



(107) 



which gives the following contribution 



G 



(2) 



\-KnmCinCiic 

in ~r Cl) 



(c?„ + c?) 



(u + cO(f + c„) (^ 



+ Ci V + Ci, 



'2CinCl 1 {Cin - Cl)^ 1 " 

{V - Cl){v - C,n) - {Cin+Cl){v-Cl){v + Cin){^-^y\ ' 

The last term, the uv correlator, has a peak (divergence) for a; < 0, a;' > 0. This feature is present in the numerical 
sinrulations of time-dependent subsonic configurations performed (together with the formation of acoustic black holes) 
in [12]. As already remarked, to have finite results we would need more realistic configurations with smooth transitions 
of width at, <Jx instead of the sudden jumps considered in our analysis. A way to take into account this is to introduce 
a cutoff in our calculation. Actually in this model, due to the two types of discontinuities we have two such cutoffs, 
J- and 
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V. HAWKING RADIATION AND DISCUSSION 



It is interesting, at this point, to ask whether step-hke discontinuities are suitable to study the Hawking effect in 
correlations and to which extent. 

In gravitational physics, the Hawking effect [l7| arises from vacuum fluctuations in the dynamical background of a 
collapsing star forming a black hole. The universality of the emitted thermal radiation (with temperature Th = ^1 , 
where k is the surface gravity at the horizon) allows us to use simple analytical models (for instance shock waves) to 
mimic the time dependence in a real process (see for instance [iMl)- The mechanism by which the black hole emits 
particles is understood in terms of pair creation in the outgoing (u) sector in the near-horizon region, one member 
outside the horizon (Hawking quanta) and the other inside (the partner, trapped inside the black hole) [l^, [20| . 
This produces a characteristic signal, a finite stationary peak in correlators (for instance, in the point-split stress 
tensor) proportional to for points situated on either side with respect to the horizon and in correspondence with 
the positions of the created quanta. 

In our context we would need a simple analytical model allowing us to study the dynamical formation of a fluid 
with a region of supersonic flow out of an everywhere subsonic configuration. As already commented in section 
4.1.1, step-like discontinuities in x can only be defined, in the hydrodynamic limit, for subsonic configurations. The 
formation of an acoustic black hole with a regular surface gravity is achieved by introducing a small transition region 
of width AcTj; around a; = where the sound velocity varies smoothly from = Ci„ (right region) to c; < \v\ (left 
region) , see Figure [51 This would give a surface gravity k ^ . 

The physics of the Hawking effect tells us to look for the uu correlations between the pairs created on the horizon 
at its formation, its members following, respectively, the trajectories x' = {cr + v)t' > (Hawking quanta) and 
X = {ci + v)t < (partner). At time t — t' this corresponds to points x, x' such that — = 0. 

A correlation of this type is present in the evolution of the 'in' left moving modes from the right, see Figure [71 j25j . 
The amplitudes A^, Aj^, depend on the details of the modes propagation in the transition region. We mention that 
an expression for them can be found by considering step-like discontinuities in the exact theory (with 5 0) and 

taking the dispersionless limit ^ — > at the end of the calculations: = ^^i:^, ~ — , = 



Ci Cr — V ^ V y Cl Cr—V 

We can then perform an analysis similar to that of section IV. C to consider the dynamical formation of an acoustic 
black hole-like configuration and taking it into account the presence of a finite Aca; by imposing a cut-off ujmax ~ ^ 
in the integration over frequencies. It is not difficult to show that the G^^-' correlator has indeed a left-right uu 
Hawking-like peak located at — and proportional to 

(2)„« ^ — n A A^ ^^^g^ 

Aumn^CrCi [v + Cl)[v + Cr) 

that has a structure similar to that found in [iJl • The next step is to extend the analysis of this paper by considering 
dispersion effects (^ ^ 0) as well. This is dealt with in detail in [Tsj . 
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HORIZON 




FIG. 6: Formation of an acoustic black hole. 
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FIG. 7: 'in' left moving modes in an acoustic black hole. 
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